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Abstract. Compton scattering of low- frequency radia- 
tion by an isotropic distribution of (i) mildly and (ii) ul- 
tra relativistic electrons is considered. It is shown that the 
ensemble-averaged differential cross-section in this case is 
noticeably different from the Rayleigh phase function. The 
scattering by an ensemble of ultra-relativistic electrons 
obeys the law p = 1— coscv, where a is the scattering angle; 
hence photons are preferentially scattered backwards. This 
contrasts the forward scattering behaviour in the Klein- 
Nishina regime. Analytical formulae describing first-order 
Klein-Nishina and finite-electron-energy corrections to the 
simple relation above are given for various energy distribu- 
tions of electrons: monoenergetic, relativistic-Maxwellian, 
and power-law. A similar formula is also given for the 
mildly relativistic (with respect to the photon energy and 
electron temperature) corrections to the Rayleigh angular 
function. One of manifestations of the phenomenon under 
consideration is that hot plasma is more reflective with 
respect to external low-frequency radiation than cold one, 
which is important, in particular, for the photon exchange 
between cold accretion disks and hot atmospheres (coro- 
nae or ADAF flows) in the vicinity of relativistic compact 
objects; and for compact radiosources. 
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Fig. 1. Angular function (in polar coordinates) for the 
Compton scattering of photons of energy E by an en- 
semble of relativistic Maxwellian electrons with temper- 
ature T e . Solid lines show Monte-Carlo simulation results. 
Dashed lines represent approximations (when available) 
given either by Eq. ( |i~6| ) or by Eq. ([H]). The dotted lines 
represent two extreme cases, one of vanishing tempera- 
ture, which corresponds to the Rayleigh angular function, 
and the other of ultra-relativistic electrons (kT e 3> m e c 2 ). 
In the latter case (presented only in the left panel), the 
scattering obeys the law p = 1 — cos a. 



1. Introduction 

The problem about the cross-section for scattering of ra- 
diation by an individual free electron, i.e. Compton scat- 
tering, is classical and its description can be found in any 
textbook on the interaction of matter and light (see, e.g., 
Bcrestctskii ct al. 1982). The Compton scattering of low- 



frequency radiation in a dense plasma for which collective 
effects are important is also well studied. In particular, for- 
mulae that describe the differential cross-section averaged 
over an ensemble of plasma-bound electrons, or, stated dif- 



ferently, the scattering angular function, are well-known 
(see, e.g., Bckefi 1966; Zheleznyakov 1996). 
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In this Letter, we would like to point out the fact, 
which is seemingly un-noticed in the literature, that the 
angular function for Compton scattering can be differ- 
ent from the classical Rayleigh function in the case of 
low-density relativistic electron gas too. We show be- 
low, providing simple analytical formulae, that the dif- 
ferential cross-section averaged over an isotropic ensem- 
ble of substantially relativistic electrons (e.g., monoener- 
getic, Maxwellian, or with a power-law energy spectrum) 
is backward-oriented, i.e. photons have a tendency to be 
scattered backwards, rather than forwards (Fig. Q). This 
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phenomenon results from the combined action of two ef- 
fects. One is that a photon is more likely to suffer a scat- 
tering from an electron that is moving toward it, rather 
than away from it (the probability of scattering is propor- 
tional to the Doppler factor 1 — cosdv/c, where v is the 
electron velocity and 9 is the angle at which the photon 
and electron encounter). The other effect is that photons 
emerge after scattering preferentially in the direction of 
the motion of the relativistic electron. The resulting angu- 
lar function contrasts the forward-oriented Klein-Nishina 
angular function which corresponds to the case of scatter- 
ing of high-energy photons on a resting electron. 

The authors certainly realize that the relativistic mod- 
ification of the differential cross-section for scattering by 
an isotropic ensemble of electrons is implicitly taken into 
account in any Monte-Carlo simulations of spectra and an- 
gular distributions of photons scattered in mildly or ultra- 
relativistic electron gases. In particular, we have verified 
that results obtained using the Monte-Carlo Comptoniza- 
tion code of Pozdnyakov et al. (1983) and the results of 
calculations that employ the formulae of this paper are in 
perfect agreement. 

Apart from a general interest, the formulae obtained 
here are easily applicable to the calculations of the scat- 
tering albedos (with respect to low-frequency radiation) 
of hot atmospheres, such as coronae or outflows above ac- 
cretion disks, ADAF accretion flows near black holes, lay- 
ers of spreading matter on the surfaces of accreting neu- 
tron stars, intergalactic gas in rich clusters of galaxies, etc. 
The issue of the modified albedo may also be important 
in studying clouds of ultra-relativistic electrons in com- 
pact radiosources, for which an external radiation may be 
an additional source of seed photons for inverse-Compton 
cooling. Our calculations of albedos demonstrate a notice- 
able reduction of the fraction of incident photons which 
can propagate inside an optically thick cloud of relativis- 
tic electrons. Therefore, less photons are capable of taking 
part in the Comptonization process, and smaller is the rate 
at which the electrons lose their energy. 

The absence of forward-backward symmetry in the 
considered regime of scattering also causes the coefficient 
of spatial diffusion of photons in hot plasma to be differ- 
ent from that in the case of non-relativistic plasma. This 
affects the shapes of Comptonization spectra and the time 
delays between soft and hard radiations coming from vari- 
able X-ray sources. A detailed discussion of both albedo 
and photon-path-lcngthcning consequences of the results 
reported here is presented in a separate paper (Sazonov & 
Sunyaev |l999al). 



2. Scattering angular function 



2.1. Scattering of low-frequency photons in 
ultra-relativistic plasma 

We start from the formula that describes the differential 
cross-section for the Compton scattering by an electron 
propagating in the direction u) with a speed v = (3c: 

i 
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(1— /3 2 ) -1 , (v, n) are the frequency and the direction 
of propagation of the incident photon, (y 1 , ft') are the 
corresponding values for the emergent photon, /i s = QQ 1 , 
and a? is the Thomson scattering cross-section. 

Consider now an isotropic distribution of electrons of 
energy 7m e c 2 (which corresponds to a given speed (3). 
We wish to average the cross-section given by Eq. ([[]) over 
this velocity distribution and thereby to find the scattering 
angular function: 



(4) 
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Here the factor 1 — /3fiw takes into account the relative ve- 
locity of the photon and electron al ong th e direction of the 
latter's motion (Berestetskii et al. 1982 ). The factor cr^ 1 
has been introduced to make p dimensionless and hence 
similar to the phase function that is used in the theory of 
radiative transfer for describing scattering processes. Ac- 
cording to the definition (||), the mean photon free path 
can be found as 

A = 1 = 1 (5) 

N e o~T J p{l, Ms) dfl'/iir N e o-T J p(7, Ms) dn s /2 ' 

where N e is the electron number density. 

We shall first consider the limiting case where 7 > 1, 
jhu/nieC 2 <C 1, i.e. the electrons are ultra-relativistic, and 
the photons are non-relativistic. The treatment can then 
be done in the Thomson limit, and, in a first approxima- 
tion, the terms proportional to hv/m e c 2 in Eqs. (^|) and 
(^|) can be ignored. The angular function will therefore be 
given by the integral 
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The main contribution to this integral is provided by 
the electrons with uj approaching Q' , which is the result 
of Doppler aberration. Therefore, we may put in Eq. (^|) 
fluj w ji a and 1 — j3£luj « 1 — ji s . The following integral 
over p! = Qu then arises: 



whereas the third term can be expanded as follows: 

7 6 (l-/VMs)(l-/W 
l + ^ 1 '^ V2 ( 1 -^ 1/2 sin. 
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The integration leads to a very simple result: 

p( 7 > l,Ms) = 1 - Ms- 
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Implementing the straightforward integrations in 
Eqs. (|l2|) and ([l3]) and keeping only the leading terms (of 
order 7~ 2 ), we derive 



The angular function (g) is presented in the left panel 
of Fig. |TJ It has a simple, apple-like shape (in polar co- 
ordinates), which is the result of the combined action of 
two effects, namely the effect of selection of electrons ac- 
cording to the incidence angle by photons and the Doppler 
aberration effect, as discussed in the Introduction. 

Using formula (^|), we can find the mean photon free 
path that corresponds to the angular function given by 
Eq. (||): A = (iV e <7y) _1 , as it should be in the Thomson 
limit. 

It is not difficult to include Klcin-Nishina corrections 
in the expression (@). We shall do that to the first order 
of accuracy. We assume as before that jhv/m e c 2 <C 1. In 
this limit, hv/^m e c 2 <C (1 — 0), and, therefore, we can 
expand the ratio v' jv given by Eq. (0) as follows: 



-l + 31n4 7 2 
H7, Ms) = 1 - Ms H j-j Ms- 
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We can finally combine Eqs. (|iTj) and (Q) to obtain 
hv 



, ... N2 -l + 31n4 7 2 
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We implemented a series of Monte-Carlo simulations to 
determine the parameter range of applicability of formula 
©:7^2, 7^^0.02m e c 2 . 

In the case of electrons obeying a relativistic 
Maxwellian distribution, which is described by the func- 



tion dN e oc 7(7 2 - l) 1 / 2 exp(-7TO e c 2 /fcT e )d7, Eq. (Jl5^ 
should be convolved with this function. The result is 
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The expression in square brackets in the equation above 
implies that corrections need to be made to the three 
terms in Eq. (^|), namely three additional terms propor- 
tional to hv/m e c 2 appear. Note that the fourth term in 
the expression X (Eq. ||) leads to a correction to the 
angular function of the order of (7/iz/ / 'm e c 2 ) 2 , which will 
be neglected here. The final result (the calculation is very 
similar to the one that had led to Eq. M) is 
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If the electron energy is not too high, say 7 <; 10, cor- 
rections of the form 7~™ to the angular function become 
important. These can be found as follows. On introducing 
spherical coordinates with the polar axis pointing along 
1~2', we obtain 



M 



Slu, = m'm s + (1 - m' 2 ) 1/2 (1 - Ms) 1/2 sine 



(11) 



Integration of the first and second bracketed terms in 
Eq. (H) then reduces to 



2(1 -/VMs) 2(1 -Ms) 



7 2 (1-/?M') 2 7 4 (1 -/3m') 3 
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where T e is the temperature of the electrons and T(a, z) = 
J x a ~ 1 e~ x dx is the incomplete Gamma function. For- 
mula ( |l6| ) is a good approximation if kT e p> 2m e c 2 , 
hvkT e <J 0.01(m e c 2 ) 2 . It is interesting to compare the first 
of these conditions with the corresponding constraint on 
the applicability of the approximation (|l5|) (see text fol- 
lowing that equation) . Evidently, the relatively poor con- 
vergence of the series (|l6|) is due to the contribution of the 
low-energy tail of the Maxwellian distribution, i.e. elec- 
trons with 7 <J (7) = 3kT e /m e c 2 . 

Various examples of the angular function correspond- 
ing to the scattering on high-temperature electrons, as 
resulted from Monte-Carlo simulations or calculated from 
Eq. (|l6|), are presented in Fig. |l|. One can see that Klein- 
Nishina corrections, which are described in the first ap- 
proximation by the second term on the right-hand side 
of Eq. Jl6|), act to reduce the probability of scattering 
in all directions. This reduction reaches a maximum at 
Ms — > — 1 (backward scattering) and monotonically dimin- 
ishes to become vanishing at p s = 1 (forward scattering). 
The effect of temperature corrections, which are, to the 
first order, described by the last term of Eq. (|l6|), is that 
as the temperature decreases, progressively more photons 
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become scattered forwards — see the pattern correspond- 
ing to kT e = 0.5m e c 2 in the left panel of Figure [l] (the 
result of a Monte-Carlo simulation is shown). 

The mean photon free path that corresponds to the 
angular function (|l6|), according to Eq. (pi), is given by 



The mean photon free path corresponding to the an- 
gular function ( |l9| ) is given by 
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This e xpres sion is equivalent to Eq. (2.17) of Pozdnyakov 
et al. (|1983| ). 

In radio galaxies, the relativistic electrons have a 
power-law energy spectrum dN e oc 7 -a o?7 with a low- 
energy cutoff 7 > j m in- The angular function in this case 
is 



wh ich is a well-known expression (see, e.g., Pozdnyakov et 
al. |1983D . 
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This formula is applicable if 7 mm 3> 1, a > 2 and 
l m inhv/m e c 2 < 1. 

12. Scattering of mildly relativistic photons in mildly 
relativistic thermal plasma (hv, kT e <C m e c 2 ) 
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The main term in this power series, 3(1 + /4)/4, is the 
usual Rayleigh function, which corresponds to the non- 
relativistic case. 

The range of applicability of formula ( |l9| ) is roughly 
kT e <j 0.05m e c 2 and hv <; 0.1m e c 2 . Some examples are 
displayed in the right panel of Fig. [l| The effect of the first- 
order temperature correction (the term oc kT e /m e c 2 in Eq. 
]l9|) on the angular function is to enhance the number of 
photons scattered at intermediate angles between 69° and 
138° (a maximum of YlkT e j (0.05m e c 2 ) per cent is reached 
at an angle of 105°) and to suppress scattering in both 
forward and backward directions (by 10fcT e /(0.05m e c 2 ) 
per cent at angles and 7r). 



